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A new Monte Carlo model, CUJET1.0, of jet quenching in nuclear collisions is applied to predict 
the jet flavor tomography spectrum in central C = — 5%, Au+Au yfs = 0.2 ATeV at RHIC 
and Pb+Pb 2.76 ATeV at LHC for fragments f = ir,D, B, e~~ fragments from quenched jet flavors 
a = g,u,c,b. The nuclear modification factors, Raa^j (pt; s,C), are predicted to exhibit a novel 
level crossing pattern over a wide transverse momentum range 5 < pr < 100 GeV that could provide 
stringent tests of jet-medium dynamics in quark gluon plasmas. 
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Introduction: Jet quenching obscrvablcs provide to- 
mographic information about the density evolution and 
jet medium dynamics in quark gluon plasmas (QGP) 
produced in high energy nuclear collisions [lj]. Exten- 
sive studies of hard probes at the Relativistic Heavy ion 
Collider (RHIC) at energies (per nucleon pair) 0.02 < 
s/s < 0.2 ATeV have recently been recently extended 
to much higher cm energies yfs = 2.76 — 5.5 ATeV at 
the Large Hadron Collider (LHC). Upgraded detectors 
at RHIC and the built in heavy quark capabilities of 
the ALICE, ATLAS, and CMS detectors at LHC will 
soon open a new chapter in jet tomography by allow- 
ing for the first time the measurement of the full jet 
parton flavor a — g,u,c,b and mass dependence of nu- 
clear modification factors, RAA-^a^f(y-,PT] s,C), for a 
wide variety of final fragments, e.g., / = n,D,B,e~, 
over much wider kinematic ranges and centrality (im- 
pact parameter) classes C. We propose in this Letter that 
the quenched jet flavor "spectrum" could provide strin- 
gent new constraints on both perturbativc QCD (pQCD) 
and string theory inspired (conformal and nonconformal) 
gravity dual holographic models of jet-medium dynamics 
in strongly interacting quark gluon plasmas. 

We report results of a new Monte Carlo pQCD to- 
mographic modelQ, CUJET1.0, that predicts a strik- 
ingly novel level crossing pattern of flavor dependent nu- 
clear modification factors at RHIC and LHC. This model 
extends the development of the GLV0, DGLV[3], and 
WHDG0 opacity series approaches by including sev- 
eral dynamical features that require extra computational 
power most easily accessible via Monte Carlo techniques. 
CUJET1.0 was developed as part of the ongoing DOE 
JET Topical Collaboration effortQ with the objective to 
develop numerical codes to reduce previous large theo- 
retical and numerical systematic uncertainties that ob- 
struct quantitative jet tomography Q and to predict new 
observables that could better discriminate between dy- 
namical models. In this letter, we focus on the jet flavor 
dependence of nuclear modification factors. 

The current CUJET1.0 code features: 
(1) dynamical jet interaction potentials that can interpo- 
late between the pure HTL dynamically screened mag- 



netic and static electric screening limits, 

(2) the ability to calculate high order opacity corrections 
up to N = 10 to interpolate numerically between N = 1 
and N = oo analytic approximations, 

(3) full jet path proper time integration over longitudi- 
nally expanding and transverse diffuse QGP geometries, 

(4) the ability to evaluate systematic theoretical uncer- 
tainties such as sensitivity to formation and decoupling 
phases of the QGP evolution, local running coupling and 
screening scale variations, and other effects out of reach 
with analytic approximations, 

(5) elastic in addition to radiative fluctuating energy loss 
distributions as well as fluctuating jet path length geome- 
tries for all flavors, 

(6) convolution over y/s and flavor dependent pQCD in- 
variant jet spectral density (without local in px spectral 
index approximations), 

(7) convolution over final fragmentation, Df/ a (x,Q), as 
well as semileptonic decay distributions. 

Motivation for the development of the Monte Carlo 
based CUJET model includes also addressing key open 
A+A phenomenology problems: (1) the heavy quar k jet 
quenching puzzle discovered at RHIC fBl llTL 121. (2) 



the surprising increase of jet transparency [13H15I ] of the 
QGP at LHC as compared to the expected linear in 
dN ch /dy(2.76, - 5%) = 16OO0 scaling of the jet opac- 
ity, as suggeste d by preliminary ALICE charged hadron 
quenching data[l6| on RpbPb-Yh+iPT] \J~s = 2.76, C = 
— 5%), and (3) the need to find hard probe observ- 
ables that can better discriminate between tomographic 
and holographic paradigms of jet medium interactions 
[ljj . As we show below the flavor dependence of LHC 
nuclear modification factor level crossing pattern out to 
Pt < 100 GeV at LHC could provide rather stringent test 
of jet- medium dynamical models. In addition, we predict 
that future RHIC jet flavor tomography out to pr < 50 
GeV could provide important consistency checks of the 
pQCD paradigm due to the large y/s variation of the 
unquenched jet distributions between RHIC and LHC. 

The CUJET Model: The CUJET1.0 code uses Monte 
Carlo techniques to compute finite order, N, in opac- 
ity contributions to the jet medium induced gluon radia- 
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tive spectrum replacing the static (Debye screened) effec- 
tive potential [Hill (GLV(eq.ll3), DGLV(cq.l7)) in the 
multiple collision gluon radiation kernel with a normal- 
ized but path dependent effective dynamical (magnetic 
enhanced) transverse momentum, q, exchange distribu- 
tion generalized from the pure HTL dynamically mag- 
netically screened modelQ to the formQ 
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where < r m = fj, m /fj, e < 1 is the ratio of a possi- 
ble nonperturbative static color magnetic screening mass, 
Mm ~ 0(g 2 T), to the standard static HTL color electric 
Debye mass, m e = gT (see also [l(|). Preliminary re- 
sults up to third order in opacity with different r m were 
reported in rcf.(BG)[l|. 

To illustrate of the proposed jet flavor tomography test 
of jet-medium dynamical models, we show here results 
only for the simplest non-trivial case with N = 1 order 
in the pure r m = HTL limitQ. At first order, the 
inhomogeneous, non-static plasma generalization of the 
dynamical MD(r m = 0) model of the induced radiated 
gluon number per (collinear) light cone momentum frac- 
tion, from massive quark jets of flavor a produced at 
position x with (transverse) energy E and propagating 
through a QGP density field, pqg_p(x, t), in azimuthal 
direction cf> is given by 
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where n a = 9C 2 (a)a 3 s /2 with C 2 (g) = 3,C 2 (g) = 4/3, 
and x+ is the fraction of plus momentum carried away 
by the radiated gluon. Also /i(t) = <?T(x(t, </>), t) is 
the local path dependent Debye screening mass with 
x(0, 0) = x as the jet production point. Here Xo.{t) = 
M 2 x 2 + +h 2 (t)(1-x + )/^/2 controls the "dead cone" effect 
due to the finite jet current quark mass for a = g, u, c, b 
as well as the local HTL asymptotic gluon thermal mass 
/.t(x + vt, t) /y/2. We include fluctuations of the radiative 
energy loss due to gluon number fluctuations computing 
P Q (e), the probability distribution of radiating a fraction 
of energy e, as in WHDG via a Poisson convolution of 
the spectrum dN g /dxE = (dx + /dxE)dN g /dx+, with xe 
the energy fraction carried away by the radiated gluon. 
We also convolute the TG model for elastic energy loss 
assuming Gaussian fluctuations as in WHDG. 

However, in contrast to WHDG, where longitudinal ex- 
pansion was only approximately taken into account via a 



mean proper time r = L/2 approximation, the CUJET 
Monte Carlo integration code computes numerically the 
full proper time integral and allows us to study system- 
atic uncertainties associated with both short time tran- 
sient nonequilibrium QGP formation effect on a proper 
time interval as well as study uncertainties associated 
with varying the decoupling times. We have tested the 
sensitivity of the nuclear modification factors to three 
extreme scenarios of the initial opacity evolution by re- 
placing the standard Bjorkcn tq/t longitudinal expansion 
factor by a temporal envelope /(t/to) that reduces to the 
Bjorken tq/t factor only after r > tq. For the present 
Letter we utilize the default form 
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where p pa rt(x)/iVp art is the normalized Glauber trans- 
verse participant nucleon density profile that depends on 
A, s, and impact parameter b, and dN c h/dy is the bulk 
charged hadron rapidity density. The level crossing pat- 
tern turns out to be surprisingly sensitive to the value of 
To as will be shown elsewhere. For illustration here, we 
set To = 0.5 fm/c. 

We assume a thermal freeze-out hypersurface T(x + 
v(0)t / ,t / ) = T f with (if = g(T f )T f = 200 McV to ter- 
minate jet path integrations. The average over initial 
transverse jet production production points is accord- 
ing to the standard binary collision Glauber Taa(x-, b = 
2 fm) density profile and the jet <j> angles are averaged 
uniformly. 

For each flavor jet and initial x, CUJET com- 
putes the fractional energy loss,e, probability distribu- 
tion P a (e;pi,x, <p) for a specified range of initial pi and 
includes delta function contributions at e = and 1. CU- 
JET then numerically averages the results over the initial 
jet production profile and initial jet <f> in order to obtain 
the final quenched partonic invariant cross section: 
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where Pi = Pf/(l — e). Note that CUJET avoids the local 
spectral index approximation R AA ~ ((1 — e Q ) 
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to 



eliminate that possible source of numerical uncertainty. 

The pQCD initial jet flavor invariant pp cross sec- 
tions, der° for y/s = 0.2, 2, 76 and y = for light 
gluon and quark (a=g,u) jet flavors were computed from 
the LO pQCD CTEQ5 code of X.N. Wang0 as in 
WHDG 11 14j. NLO charm and bottom quark invariant 
cross sections for both RHIC and LHC were provided by 
R. Vogt 0. 

RHIC and LHC Results: 

Our central physical tomographic assumption leading 
to our main result, Figure 1, is that aside from the 
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FIG. 1: Illustration of jet flavor tomography level crossing pattern of nuclear modification factors versus pt at y = for 
ir, D, B, e fragmentation from quenched g, u, c, b jets in Au+Au 5% at RHIC (left side) and extrapolated to Pb+Pb 5% at LHC 
(right side) computed with the dynamic CUJET1.0 model at leading N = 1 order in opacity. The opacity is constrained at 
RHIC by a fit to a reference point, R AuAu (p T = 10 GeV) = 0.2, with dN/dy(RHIC) = 1000 and a a = 0.3. The extrapolation 
to LHC assumes dN c h/drj scaling of the opacity as measured by ALICEflJ. The level crossing points are denoted by circled X. 
Note the possible inversion of n, D, B levels predicted by CUJET at high pt at LHC and a partial inversion at RHIC arising 
from competing dependences on the jet parton mass of energy loss and initial pQCD spectral shapes. The crossings points and 
level spacings are sensitive to dynamical assumptions as will be reported elsewhere. 



unavoidable dependence of the initial pQCD par- 
tonic invariant cross sections, the sole yfs dependent nu- 
clear input in CUJET is the variation of the bulk fi- 
nal pion rapidity density, dNAA/dy with t/s. There- 
fore, as in WHDGllflil]. we assume that the charged 
particle pseudo-rapidity density dN^ c /drj = 1600 re- 
ported by ALICE [l3 for C = - 5% central Pb + Pb at 
2.76 ATeV translates directly into a factor 2.2 increase 
of pqgp a t LHC relative to RHIC (for the same cen- 
trality). Furthermore, we fix the RHIC partonic level, 
dN/dy = 1000 with a s = 0.3, in Eq.(3) to constrain one 
reference = 10 GeV point of pion R\a = 0-2 at RHIC. 
The RHIC constrained extrapolation to LHC is then pa- 
rameter free (assuming a s , tq and Tf do not vary with 
y/s). Systematic uncertainties associated with variations 
of the effective jet medium coupling 0, [H| as well as ini- 
tial and final temperature field profiles will be reported 
elsewhere. Because of CUJET includes dynamical MD 
magnetic enhanced potential [8|, our assumed moderately 
small coupling is sufficient to account for the RHIC data 
with dN(RHIC)/dy = 1000 as we check below (see fig. 
2). 

In Fig.l on left side, the splitting between pion and 
electron Raa is found to remain large below 10 GeV 
in spite of the use r m dynamically enhanced potential. 
Also shown are the separate contribution of D and B 
mesons that future detector upgrades at RHIC could 
test. The importance of experimentally isolating and 



observing charged heavy mesons cannot be overstated 
since the mass splitting between c and b jets is a partic- 
ularly robust prediction of pQCD in a deconfincd QGP 
medium. The heavy quark quenching puzzle therefore re- 
mains unsolved in the CUJET calculations shown. Holo- 
graphic gravity dual models predict a qualitatively dif- 
ferent heavy quark quenching pattern [18j |. The novel 
inversion of the 7r < D < e < B Raa hierarchy ordering 
at high pt is due to the steeper initial invariant jet dis- 
tributions of c and b jets at RHIC[20T] as also noted in 
21 1 . 

Future LHC data will provide an opportunity to test 
details of the jet flavor dynamics at densities at least a 
factor of two greater than at RHIC. On the right side of 
Figure 1, the RHIC constrained LHC extrapolated jet fla- 
vor tomographic pattern is shown. With the much wider 
kinematic window accessible at LHC, the predicted flavor 
dependent pt spectrum of nuclear modifications is seen 
to involve multiple level crossings that are qualitatively 
different than at RHIC energies because of the compli- 
cated interplay between flavor dependent spectral shapes 
and opacity enhanced jet energy loss. While the details of 
the level crossings depend sensitively on specific dynami- 
cal details such as the effective coupling, initial formation 
time, and freeze-out, that the differential in pT,y/s,C, 
and A variations of the jet flavor quenching pattern ap- 
pear particularly promising to discriminate between jet 
medium dynamical models already at this single inclu- 
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FIG. 2: Dynamically (r m = 0) enhanced jet quenching[8| at 
RHIC (blue dashed) and LHC (red solid) with the CUJET1.0 
at pt = 10 GeV RHIC as in Fig.l. Central - 5% PHENIX 
7T° RHIC data[!2] (blue dots) and preliminary ALICE h± 
LHC data[l(| (black dots within pink systematic error band 
that take into account preliminary pp 2.76 reference data) are 
compared to predictions. For comparison, previous LHC pre- 
dictions from (IV02)[l| (grey) and WHDg[I3 (purple) based 
on static Debye screened jet medium interactions are shown 
as well. 



sive level. Of course flavor tagged dihadron and hadron- 
lepton correlations will enable much more detailed quan- 
titative information to be extracted eventually. 

In Fig. 2 we compare our CUJET constrained 

predictions for pion R\ a (pt) to central 0-5% 

PHENIX/RHICHI and preliminary ALICE/LHC0 
data. We also compare to previous LHC predictions 
that used the static Debye screened potential (IV02)[l| 
and recently updated WHDGllQ models. CUJET is 
seen interpolate between the later two curves, and the 
curves have a tendency to fall below the preliminary 



LHC data (as do several other recent predictions [13lll5| 
not shown). This suggests the intriguing possibility that 
the jet-medium coupling at LHC is possibly weaker than 
at RHIC. See (T^Ji however, for a detailed discussion of 
many open caveats in this connection. It also remains 
an important open question if bulk elliptic QGP flow at 
LHC shows any signs of weaker coupling [23T], 
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